Parton Transport via Transverse and Longitudinal Scattering in Dense Media 
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The effect of multiple scatterings on the propagation of hard partonic jets in a dense nuclear 
medium is studied in the framework of deep-inelastic scattering (DIS) off a large nucleus. Power 
counting arguments based on the Glauber improved Soft-Collinear-Effective-Theory are used to 
identify the class of leading power corrections to the process of a single parton traversing the extended 
medium without emission. It turns out that the effect of longitudinal drag and diffusion (often 
referred to as straggling) is as important as transverse scattering, when relying solely on power 
counting arguments. With the inclusion of momentum exchanges in both transverse and longitudinal 
directions between the traversing hard parton and the constituents of the medium, we derive a 
differential equation for the time (or distance) evolution of the hard parton momentum distribution. 
Keeping up to the second order in a momentum gradient expansion, this equation describes in- 
medium evolution of hard jets which experience longitudinal drag and diffusion plus the transverse 
broadening caused by multiple scatterings from the medium. 



I. INTRODUCTION 

The modification of hard jets in hot (and or dense) 
extended media is a forefront topic in relativistic heavy- 
ion collisions as it may provide a sensitive measure of 
the properties of the produced highly excited matter 0- 
Q . The partonic showers initiated by these hard off-shell 
patrons are modified due to the scattering experienced 
by each of the partons in the course of the propagation 
of the cascade through the medium. Such a modification 
is confirmed by multiple experimental observations such 
as the significant depletion of high transverse momen- 
tum (pt) hadron production compared to that in binary- 
scaled proton-proton collisions at the same energies p-Q ■ 

The cause of this modification is a combination of drag, 
diffusion and, most importantly, stimulated emission ex- 
perienced by the shower while propagating through the 
medium |9l-[l9|. Sophisticated quantitative jet quench- 
ing calculations have been performed for single inclusive 
hadron suppression [2(il - l25j . as well as di-hadron corre- 
lations [26l - [29| and photon- hadron correlations j30l - l32j . 
Following a successful description of the data, a signifi- 
cant amount of effort has now been placed on the quan- 
titative extraction of jet transport parameters, such as 
q = d(ApT) 2 /dt and e = dE/dt, in hot and dense media. 

The current manuscript represents an extension of 
the theory of parton energy loss, based on the Glauber 
improved Soft-Collinear-Effective-Theory (GiSCET) ap- 
proach [33|, Hi| . This approach is closely linked with the 
multiple scattering versions of the Higher Twist (msHT) 
scheme [l9l l35j and our presentation will lie in-between 
these two approaches: we will not use the GiSCET La- 
grangian; however, the power counting rules of the differ- 
ent momentum components will be similar to those used 
in Ref. H3- 

There is one exception to the above statement which 
constitutes the raison d'etre of this paper. In all pre- 
vious attempts to compute the effect of a medium on 
the hard jet parton, the medium gluons are assumed to 



have vanishingly small light cone components, both in 
the direction of the jet's momentum and in the direc- 
tion opposite to it. In this work, as in Refs. [l9|, l35j . 
the struck parton is assumed to possess a momentum 
q = (q + ,q~ ,q±) ~ (A 2 Q, Q, XQ), where Q is a hard scale 
and A is a small dimensionless number. In the previous 
papers, the gluons off which the jet scatters were assumed 
to have a momentum k ~ (A 2 ,A 2 ,A)Q- The transverse 
momentum components are required to have the same 
scale as those of the jet. The (+)-componcnts have to 
be ~ A 2 Q to maintain the virtuality of the jet. How- 
ever, there is no physical reason why the (— )-component 
of k has to be ~ A 2 Q. Indeed it may be as large as XQ 
without changing the virtuality of either the hard par- 
ton or the Glauber gluon (there will be a difference in 
the effect on the medium). In this paper, we explore the 
fate of a hard parton with propagates through a dense 
extended medium interacting with such transverse and 
longitudinal Glauber gluons. 

In the current effort, the sole focus will lie on the effect 
of multiple scatterings; radiation and radiative energy 
loss for the propagating hard parton are not included in 
this work. The longitudinal momentum exchange during 
multiple scatterings leads to the drag and diffusion of the 
propagating jet, which has often been understood as the 
collisional (or elastic) jet energy loss. While the relative 
importance of the two mechanisms is still in dispute, col- 
lisional energy loss cannot simply be neglected. This is 
true, not only for the suppression of single inclusive light 
and heavy hadrons [2(| Hj], but also for jet shower evolu- 
tion, energy loss distribution within and outside the jet 
cone, as well as energy and momentum deposition into 
the medium by the jet shower [36l- 38] . The longitudinal 
momentum exchange may affect induced radiation from 
the hard parton as well; this will be presented in an up- 
coming publication. The paper is organized as follows. 
In the next section, we present the three dimensional 
momentum distribution of the outgoing hard parton at 
leading twist. In Sec. Ill, we investigate the effect of mul- 
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tiplc scatterings on this distribution. A class of the higher 
twist diagrams which are length enhanced arc identified. 
Then the gradient expansion is invoked and the number 
of multiple scatterings is resummed. We finally derive a 
differential equation for the time evolution of parton mo- 
mentum distribution as affected by multiple scatterings. 
Sec. IV contains our summary. 



where the J^x runs over an possible final states except 
the stuck hard jet. J M = Qqipq^ipq is the hadron elec- 
tromagnetic current, with Q q the charge of a quark of 
flavor q in units of the positron charge e. Here our focus 
is on the final state interaction between the medium and 
the stuck quark and thus the discussion will be centered 
around the hadronic tensor; the leptonic tensor will not 
be discussed further. 



II. LEADING TWIST AND PARTON 
DISTRIBUTION FUNCTIONS 

Consider the process of semi-inclusive DIS off a large 
nucleus where one hard parton with a momentum l q is 
produced, 



e(L 1 ) + A(P A )^e(L 2 )+q(l q )+X 



(1) 



In the above expression, L\ and L 2 represent the mo- 
menta of the incoming and outgoing leptons. The in- 
coming nucleus of atomic number A has a momentum 
Pa = Ap, with each nucleus carrying momentum p. In 
the final state, all hadrons are detected and their mo- 
menta are summed to obtain the hard parton's momen- 
tum l q . X denotes that such process is semi- inclusive. 

Throughout this work, we utilize the light-cone com- 
ponent notations for four vectors (p = [p + ,p~ ,p±\) with 

p+ = (E+ Pz )/V2,p- = (E- Pz )/V2. (2) 

In the Briet frame of DIS, the incoming virtual photon 
7* has a momentum q and the nucleus has a momentum 
Pa, 

q = L 2 -L 1 = [-x B p + ,q~ ,0 ± ], P A = A[p+ ,0,0], (3) 

where xb — Q 2 /(2p + q~) is the Bjorken variable in this 
frame. 

The double differential cross section of the semi- 
inclusive process in which a jet with a momentum l q is 
produced may be expressed as 



El 2 cI<j 
<PL 2 dH q 



a e 1 T 
2irsQ± ^ dH 



(4) 



where ,s = (p + L\) 2 is the total invariant mass of the 
lepton-nucleon collision system. The leptonic tensor is 
give by 



The semi-inclusive hadronic tensor is defined as 



(5) 
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FIG. 1: Leading twist contribution to the semi-inclusive DIS 
hadronic tensor 

The leading twist contribution to semi-inclusive DIS 
is obtained by expanding the products of the currents 
at leading order in a s , as shown in Fig. [1] This rep- 
resents the process where a hard quark, produced from 
one nucleon of the nucleus, exits the nucleus without fur- 
ther interaction. The leading twist contribution to the 
semi-inclusive hadronic tensor may be expressed as 



x B ). 



(T) 



In the above expression, The projection tensor gj" is de- 



fined as g^ = g» v - g»~ g u+ - g^ g v ~ . The factor C A 
represents the probability to find a nucleon state with 
momentum p inside a nucleus with A nucleons. The func- 
tion f q (x) represents the parton distribution function of 
a quark with flavor q in a nucleon with a fraction x of 
the forward momentum p + of the nucleon, 



/«(*) 



2ir 



-ixp y a 



(myo^mip), (8) 



where j + and 7 are used to project out the leading 
hclicity components, 
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(9) 



x 



{2n)' l S 4 (q + P A -Px) 



x(A\j^o)\x)(x\r(o)\A). 



(0) 



In the above expression, \A) represents the initial state 
of an incoming nucleus with A nucleons with momen- 
tum p per nucleon, averaged over spins. The state \X) 
represents the general final hadronic or partonic state, 



In the limit of very high forward momentum, the in- 
coming parton carries vanishingly small transverse mo- 
mentum (pj_ ~ XQ), the differential hadronic tensor for 
the momentum distribution of the final quark may be 
approximated as, 



dW, 



A^LU 



d 2 l q ±dl q 



= Wr 



Anvr,,- 



5(1- - q-W{l qA _). 



(10) 
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In the following section, wc will investigate how the mul- 
tiple scatterings from the medium change the momentum 
distribution of the final quark. 

III. MULTIPLE SCATTERINGS AND PARTON 
EVOLUTION IN MEDIUM 

In this section, the effect of multiple scatterings on the 
3-dimcnsional momentum distribution of the final quark 
is presented. The higher twist contributions are obtained 
from those diagrams which include the expectation val- 
ues of more partonic operators in the medium. Such con- 
tributions are usually suppressed by powers of the hard 
scale Q 2 , but a sub-class of these contributions may be 
enhanced in extended media due to the longitudinal ex- 
tent traveled by the struck quark. 

The generic diagram being considered here is shown 
in Fig. [5] which describes the process that a hard virtual 
photon strikes a hard quark in the nucleus with momen- 
tum p' (po in complex conjugate) at location j/q = 
(yo in complex conjugate). The stuck quark is then 
sent back through the nucleus and has momentum q[ 
(qi in the complex conjugate). During its propagation 



through the nucleus, the hard parton scatters off the 
gluon field within the nuclear medium at locations y'j 
with < j < m (yi in the complex conjugate with 
< i < n). In this effort, we are only considering the case 
where the hard quark propagates through the medium 
without radiating. Radiation will be dealt with in a fu- 
ture effort. Through each scattering, the hard parton 
picks up momenta p'j (pi in the complex conjugate). Us- 
ing the momentum conservation at each vertex, we may 
denote the various momenta in the picture as follows: 



qi+i = q% + Pi = q + ^Pj = q + Ki, 

3=0 

i 

q' l+ i =ql+p' i = q + Y,p' j = q + K (n) 

3=0 

where for convenience we have defined new momentum 
variables Ki = J2]=oPj' an d %l = Y?j=oPj> which rep- 
resent the total accumulated momentum exchanged be- 
tween the hard quark and the nuclear medium. For such 
a diagram, we may write down the hadronic tensor as, 



w Ay - 



ran / . 
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y q g N 



n 



d 4 qi 



-iqi\vi-\-Vi) 



-u-(v-n-y' m ) 




Wj-iv'j-y'j-i) 



(AmyoW I 



i -qi 




7 ^(0)| A). 



(12) 



Here, Q q is the electromagnetic charge of the quark in units of an electron's charge, N c is the number of colors, 
and T ai (T a i ) represent Gcll-Mann matrices. To simplify the above expression, we may change the integral variables 
<7i + i — > pi and q'j+i — > p'j, and introduce the n th exchanged momentum p n in the complex conjugate by inserting the 
following momentum conserving ^-function, 



d A p n 
{2-kY 



(27T) 4 6 4 (l q ~q + K n ). 



(13) 



3=0 



The m th momentum p' m in the amplitude is determined by all other momenta, i.e., p' rn = 
The phase factor appearing in the above expression can be collected as (n"=o e ~ Wi Vi ) (jl 

Before progressing further, we outline the power counting relations which underly the various approximations 
that will be carried out in the remainder of this paper. The momenta exchanged with the medium has transverse 
components that scale as k± ~ XQ where Q represents the hard scale in the problem. Different from prior calculations 
we will also insist that k~ ~ XQ. Given the large Lorentz boost, the (+)-component of k may be quite large, however, 
due to the requirement that the propagating parton not go off-shell by more than XQ, we will obtain that k + ~ X 2 Q. 
We now consider the propagators after each scattering. For the quark momentum after i th scattering, we have 



q}+i 



{q + Ki) 2 = 2p + q {l + Kjq )[-X B +X i -X Di 



(14) 
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FIG. 2: An order of n + m contribution to hadronic tensor W^ v with n gluon insertions in the complex conjugate and m gluon 
insertions in the amplitude. 



where we have defined a few momentum fraction variables, 



Xl = VlEj = V -1- = X Di = Y] X Di = — T- ~ — . -• (15) 

j^Q J^qP P J^o 2p+q-{l + K l /q~) 

We may collect the contributions from all the internal quark line denominators together with the on-shell delta 
function for the final quark line l q as, 



„ in-i ^ \ im-i 



where the factor C q stands for, 

nirw)(nrnkF)- (17) 



U=0 i / " / \ j=Q 



In so doing, we have retained the first correction in A stemming from the longitudinal loss experienced from exchanges 
with k~ ~ XQ. As argued in Ref. [Hj], in the high energy (Q —t oo) and collinear (A —> 0) limit, we may drop the 
transverse parts of the 7 matrices which connect the quark current to the virtual photon, which implies, 

(4>(y)dm) « 1 —{i>{y) 1 - Y 6m)- (is) 

Also, noting that A + ~ A 2 Q and A± ~ A 3 Q, in A - = gauge we approximate, 

1 -A{y)K 1 -A + {y). (19) 

One should point out, once again, that in this manuscript we will demonstrate that the effect of longitudinal drag 
is the same order as that of transverse scattering; however, in the derivation, we will retain corrections up to order 
A. While the order A corrections to the propagators are retained in Eq. (fT7|) . we neglect contributions to the vertices 
from A± in Eq. (|19[) . This is due to the fact that the first non-zero contribution containing an Aj_ term is of the from 
7-1 ' 9i±7-L ' Aj_ which is of order A 2 . 
With these simplifications, we obtain 

(A\<P(y )r ( f[ 7 . qci ■ A^{ Vi ) j 7 ■ U J} 7 • ^(y-)7 ■ ] 7^(0) I A) 

\i=l / \j=m J 



{WiVo)^ II^ +a '^) II A+ai tti) ] ^(0)|A)Tr 

\i=l / \7=m 



J_ 

2 



7^ ( n 7 ■ ) 7 ■ i i n 7 7 ■ ^ 1 7 ' 

\i— 1 / \ j— m 



(20) 
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We now simplify the trace over spinor indices. Using qi+i = q + Ki = q + X^j=oP*' ^ ne trace reads, 
f[ 7 (q + IU-lhA J-(q + K n ) [ H 7-7 • (q + K'^) 



T = Tr 



2 ' 



\i=l 



(21) 



where 7 • (q + Ki) = j + (q~ + K^) + -f~(—XBP + + K^) — y± ■ Ki t ±. The above trace can be easily carried out by 
noting that 7~7 _ = 0, as a result, all terms in 7 ■ (q + Ki) containing a 7" vanish. Terms containing 7j_ • K± can 
only be included along with a j±A± and thus contribute A 2 corrections to the cross section. As a result, we obtain, 



T=H?f) 



_ \ n+m+1 



(22) 



Substituting the above simplifications into the hadronic tensor, one obtains, 



\3=m 



d% 



K n ) 



n/^)(nH(n/l/^)(ff 



dx'j 
"2V 



V,i=0 



(27r)(5(-2; B + x n - x Dn ) ] 



m — 1 



JJ e ix 'tp + (y'j -y' m ) e i P , j-(y'j-y' m ) 



-x B + Xi - x D i - «e 




1 



3=0 



-Xg + X 1 j — x' D . + it 



(- g ^)(A\iiy )l- n A+a, (^) ] 



V>(o)|A), 



(23) 



where we have performed the integration over / + using the delta function (2tt)6(1 + — q + — K^), and have changed 
integration variables pf — > Xi = p^ /p + , p^ — > x'j = /p + . We have introduced a vector notation for momentum 
and space coordinates for convenience: p = (p~,p±) and y = (y + ,y±), and their dot product represents ft - y = 
p~y+ -p ± ■ yj_. 

Now we perform the integration over the momentum fractions Xi and x'a- The integration over x n is performed 
using the last delta function which forces the cut line to be on shell and yields, 



+ x B + xd, 



i=l 



XB + XDr 



The (+)-component of the phase factor now reads, 

/n—l \ 



m—1 



3 -t(i fl +i D „)p i + y n 



n 



n^ +b7 " y ' ;) | =r+r+. 

.0=0 



(24) 



(25) 



In the above equation, and represent the two phase factors associated with the Xi integration and the re- 
integration, respectively. The remaining integrations over the momentum fractions Xi and x'j may now be performed 
using contour integration. One starts from the propagators adjacent to the cut and proceeds to the initial hard 
scattering vertex. In the following, the integrations over the fractions in the complex conjugate (x^s) will be described 
in detail; the integrations over the fractions Xj's in the amplitude are completely analogous. 

The first integration is the one over the momentum fraction x n —i. Isolating the piece related to x n -± integral, we 
can perform the integral by closing the contour of with a counter-clockwise semi-circle in the upper half of x n -\ 
complex plane and obtain 



dx r 



-ix n - 1 p+(y n l -y n ) 



2n -X B + X n -l + X n -2 - XDn-1 ~ « e 



,,- \p-i(-Xr l -2+XB+XD n -l)p + (y n _ 1 -V n ) 



(26) 
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Here the ^-function represents the physical effect that the quark line is propagating from y n _i to y n . Correspondingly, 
the phase factor combined with the above from the contour integration becomes, 



iP (Vi -Vn-l) 



(27) 



The phase factor is found to have a structure similar to the case without longitudinal drag [35[. Thus, the remaining 
integrals over the longitudinal momentum fractions Xj's can be done in a similar way. Finally, we obtain 



n« 

\i=i 



-ix Di p y i 



Il%r-%-i) 



(28) 



\i=l 



The integrations over the momentum factions x'j's in the amplitude are completely analogous, except that a factor of 
(— i) instead of i is associated with the 9 function. Such difference originates from the contour integral of x n -\ with 
a clockwise semi-circle in the lower half of the complex plane. 

After having performed the integrations over the internal quark lines, we may combine various parts and obtain 
the hadronic tensor as follows, 




I J J ^Pj-Vj 



3 = 1 



(29) 



The expression derived above is quite general in the sense that we have not made any assumption regarding the nature 
of the nuclear states. In the following, we will make specific assumptions on the nuclear states to simplify the above 
expression. Here we only consider the case of n = m, i.e., the same number of scatterings in both the amplitude and 
the complex conjugate. The case where n ^ m, either constitute higher twist nucleon matrix elements or unitarity 



corrections to terms where the produced quark experiences scatterings min(ri, m) times [35 



For the case of 2n (n = m) gluon insertions, we may first simplify the matrix elements of the gluon vector potentials 
in the nuclear state. In this work, the nucleus is approximated as a weakly interacting homogenous gas of nucleons. 
Such approximation is sensible at very high energy, where nucleons appear to travel in straight lines and are almost 
independent of each other over the time interval of the interaction of the hard probe due to time dilation. As a result, 
we may decompose the expectation of field operators in the nuclear states into the expectation in the nucleon states. 
As a nucleon is a color singlet, any combination of quark or gluon field strength insertions must be restricted to a 
color singlet. Therefore, the first non-zero (and the largest) contribution comes from the terms where 2n gluons are 
divided into singlet pairs in separate nucleon states where one gluon of the pair is in the amplitude and one is in the 
complex conjugate. 

Using the time ordered products of 9 functions above, the expectation of the field operators in the nuclear state 
may be decomposed as, 



IP 




(pmy )^m\p) 



Y[(p\A^(yAA+<(y>)\p) 



(30) 



In the above decomposition, we keep only the largest contribution arising from the terms where the expectation of each 
pair of partonic operators is evaluated in separate nucleon states. The yi integrations arc carried out over the nuclear 



volume under the constraints imposed by the string of 0-functions. The factor C p 



represents the probability 
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to find n + 1 nuclcons in the vicinity of the positions yo, JJx • • ■ y n - In the case of non- interacting nucleons in a nucleus 
with a uniform density, this factor may be approximated as 

where p is the nucleon density inside the nucleus, and the factor l/(2p + ) originates in the normalization of nucleon 
state. Note that in the case of a nucleus which docs not have a uniform density the coefficient will possess a spatial 
dependence. One may average over the colors of the gluon fields, 

(A a (y)A»(0)) = J ^- J (A(y)A(0)). (32) 

The average over the colors of quark field has already been carried out to give the factor 1/N C . Now the overall trace 
over the color factors reduces to 



1 1 -Tr 



J T ai \ JJ T CH 

i=l / \j=n 



N*-l 



(33) 



N c (iV c 2 - 1 
We may now obtain the hadronic tensor as 

J dy - J d 3 y Q J ^e-«>-*e- fa ^ + ''o(p|^( I/b )^V(0)IP> 

(ft / d »r / <V* 7 %r - Vi-jetf: wT-i)J 

(ft/^W / ^i/ |^ / ^ e -®*e^^e-^^ + ^e-°^ + ^(p|A+(y,)A+(^)b)j , (34) 

where we have changed variables from dpg to dp' n . 

The above expression may be further simplified by changing variables (yi,y'j) — > (Yi,Syi), 

Yi = (vi + Vi)/2, Sy i =y i -y' i . (35) 

For a large enough nucleus one may insist on translational invariance for the expectation values of gluon operators, 
i.e., 

(p\A + ( yi )A + (yl)\p) ~ ( P \A+( yi - yl)A+(0)\p) ~ (p\A + (6 yi )A + (0)\p). (36) 

Now only the phase factor depends on the average values Y+, thus we may perform the integration over the phase 
factor and obtain 

d 3 y z I d^e-^-tf'e*^ = / d 3 Y t [ ^fae-*®*- fo** e -W* = (27r) 3 <J 3 (& / cPfae-*****. (37) 



The above delta function will set the momentum fractions = a;^-, .. For the product of ^-functions, since both Sy i 
and Sy~_ 1 are within the nucleon size (small compared to the nucleus size Y~), we may simplify it as, 

%r - vT-iWi - v'7-i) = 9{Yr - Yrj. (38) 
The time-ordered product of ^-functions gives, 



(39) 



where L is the extent of the nucleus size. Combining the above simplifications and using the delta function to 
perform the integration over p' iy wc obtain the differential hadronic tensor as, 



dW At±v 



(p\A+(6 yi )A + (0)\p) 



d 3 i q ,t: V Jo j j j w jv c - 1 ) 

( II e- teD ^ +i »«" ) 5 3 ( f, - g - f> ) . (40) 
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With all the simplifications carried out above, it is now possible to carry out the resummation over multiple 
scatterings. As in the case without longitudinal exchange, we will not resum the entire hadronic tensor as outlined 
above. Instead we introduce further simplifications arising from the collinear approximation, which restricts the 
exchanged momenta to be small compared to the energy of the hard parton. One Taylor expands the hard part (the 
last line in the above expression) in terms of the exchanged momenta pi around pi = 0, 



*= n 



1 



dpi 



1 a fi _d_±_ 

2 PiPi dpf dp P 



(41) 



In the above equation, the values of a, (3 take the values of "— " and "_L" . Here we only keep the terms in the 
expansion up to the second order derivatives, with the assumption of small momentum exchange in each of the 
multiple scatterings. The extension to include higher derivative terms should be straightforward, which corresponds 
to including higher moments of the momentum distribution of the exchanged gluon. In the above expansion, these 
terms involving the zeroth order term (without derivative) represent the case where the exchanged gluon carries zero 
momentum. Such terms corresponds to gauge corrections to diagrams with lower number of scatterings and will not 
be considered further. These terms should be understood to be included in the gauge invariant definition of the 
transport coefficients. 

Using integration by parts, the factors of exchanged momentum may be converted into the appropriate derivatives 
over relative position, 



-m-S& {plA+ (Q) lp)p a _d_ = e _ ift . m ( _ . } {pl g a A+ (Q) ^ _d_ ^ 



-^•^(p|^+(^)A+(0)b)pfpf 



d P f dp p 



-^■ s v*(p\d a A + (5ui)d p A + (0)\p) 



_d d_ 

dpfdpj' 



(42) 



Since the Taylor expansion imposes the condition pi — > on the hard part H, it now no longer has functional 
dependence on pi. Therefore, we may perform the integrations over pi and Siji. In this effort, we only keep the leading 
term arising from the delta functions when evaluating the derivatives on the hard part and ignore the terms from the 
derivatives on the phase factor e~ lXDiP Sy * which result in spatial moments of the two gluon matrix elements such 
as (A + (5y~, Sy±)Sy~ A + (0)). With the above simplifications, the differential hadronic tensor now reads, 



dH a 



W, 



AjjLV 



1 



n 



dYi 



n 9 ! n d2 1 
-u L1 -— + t; u l2^-= + x 

op i l v Vi * 



D T2 Vl 



,o.(43) 



In writing the above expressions, we have focused on the case with unpolarized initial and final states and only kept 
terms with non- vanishing coefficients by utilizing the symmetry of the system. The jet transport coefficients Dl\, 
Dl2 and Dt2 in the equation above are defined as, 

Dli =a 2 Jjf— l J dy-^(p\id-A+(y-)A+(0)\p), 



r> 2 C F 



D 



T2 



' N?-l 



dy~-? T (j 3 \d-A+(y-)d-A + (0)\p), 
2p + 

dy-/-{p\d^A+{y-)d^A + {0)\p). 
2p + 



(44) 



One can immediately see that if the longitudinal momentum exchange is the same order as transverse momentum 
exchange, then Dl2 ~ Dt2 as stated in the Introduction. This indicates it is important to include both effects 
when considering the evolution of hard jets in dense nuclear medium. In what follows, we will demonstrate that 
these coefficients D^i, and D?2 are, up to an overall factor, related to elastic energy loss rate e, the diffusion 
in longitudinal momenta e 2 and the diffusion in transverse momenta q. We note that the above definitions of jet 
transport coefficients are not gauge-invariant. The gauge invariant definition of the transport coefficients can be 
found with the inclusion of higher order terms where k± <C XQ. The summation over an arbitrary number of such 
soft gluon insertions leads to the appearance of Wilson lines between the the gluon field strength operators, which 
renders the operator product gauge invariant. For the case of jet broadening, more detailed discussion can be found 
in Ref. 0,E3|. 

The resummation over an arbitrary number of multiple scatterings may now be performed, 



dW A ^ 



E 



dw£r 

d 3 L 



(45) 
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where we have denned the final quark distribution function <p(L , L , l q ±) 




dY; 



-D 



Ll 



_d_ 

dp; 



i d 2 



S 3 \ 



lpi---Pn=0- 



The quark distribution function, after resumming over an arbitrary number of scatterings, reads, 

d 1 d 2 1 
- + -D L2 h -i 

2 L2 8 2 la 2 



exp L 



D 



Ll 



~rD T2 Vf 



5(l--q-)S 2 (l ql _), 



(46) 



(47) 



where we have converted the derivatives over pi to the derivatives over l q . 

It is easy to see that the above distribution function <p{L~ , l~ , l q j_) for the final quark satisfies the following evolution 
equation, 

^ + \d L2 ^ + \d T2 VI^ <P(L-,l-,l q± ). (48) 

Such an evolution equation describes the evolution of the three-dimensional momentum distribution of propagating 
partons which suffers multiple scatterings without radiation. The three terms in the above evolution equation represent 
the contributions from longitudinal momentum loss and diffusion, and the diffusion of transverse momentum. Using 
the initial condition: 4>(L~ = 0,l~,l q ±_) = 5{l q — q~)6 2 (l q ±_), the distribution <fi has the following solution, 



dL- 



D 



Ll 



,-(l--q-+D Ll L-) 2 /(2D L2 L~) -l 2 q± /(2D 7 



y/2TrD L2 L- 2nD T2 L- 
From the above solution, one may obtain 

(l~)=q- -D L1 L~, ((l-) 2 )-(l-) 2 = D L2 L- 1 (l 2 )=2D T2 L 



(49) 



(50) 

The coefficients Dli, Dl 2 and Dt 2 are defined on the light cone, and they are related to elastic energy loss rate 
e = dE/dt, the diffusion of elastic energy loss e 2 = d(AE) 2 /dt, and the transverse momentum diffusion rate q = 
d(Ap T ) 2 /dt as 



e = D L i, e 2 = D L2 /V2, q = 2V2D T2 . 



(51) 



IV. SUMMARY 

The current manuscript forms the first of a series of 
efforts meant to incorporate elastic drag and scattering 
induced emission into one seamless formalism. In this 
work, we have studied the effect of multiple scattering 
on the propagation of a single hard parton in the dense 
nuclear medium. The developed formalism will be ap- 
plied to both jet modification in DIS on a large nucleus 
as well as in the case of heavy-ion collisions. To study the 
effect of multiple scattering on a single patron, we have 
suppressed radiation from the parton. This will be dealt 
with in a future effort, which will include the scattering 
of both the parent parton as well as the radiated gluon. 

A class of higher twist corrections which are enhanced 
by the large extent of the nucleus are resummed to ob- 
tain the length dependent three dimensional momentum 
distribution of a hard quark which exits after travers- 
ing the dense medium. We include both transverse mo- 
mentum broadening as well as the drag and momentum 
diffusion in longitudinal direction between the hard jet 



and medium constituents simultaneously when consid- 
ering jet propagation in a dense medium. Momentum 
power counting based on SCET-Glauber scaling indicates 
that both these corrections are of similar size, though the 
longitudinal drag is still very small compared to the en- 
ergy of the hard quark. A differential equation has been 
derived for the time evolution of jet momentum distri- 
bution as affected by the multiple scatterings from the 
dense medium without radiation. The three coefficients 
involved Dl±,Dl2, Dt 2 , have been related to the elastic 
drag e, longitudinal straggling e 2 and transverse momen- 
tum diffusion coefficient q. The combined effect of all 
these three coefficients on the radiated gluon spectrum 
will be explored in a future effort. 
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